Kondo temperature and screening extension in a double dot system. 
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In this work we use the Slave Boson Mean Field Approximation at finite U to study the effects of spin- 
spin correlations in the transport properties of two quantum dots coupled in series to metallic leads. Different 
quantum regimes of this system are studied in a wide range of parameter space. The main aspects related to 
the interplay between the half-filling Kondo effect and the antiferromagnetic correlation between the quantum 
dots are reviewed. Slave boson results for conductance, local density of states in the quantum dots, and the 
renormalized energy parameters, are presented. As a different approach to the Kondo physics in a double dot 
system, the Kondo cloud extension inside the metallic leads is calculated and its dependence with the inter-dot 
coupling is analyzed. In addition, the cloud extension permits the calculation of the Kondo temperature of the 
double quantum dot. This result is very similar to the corresponding critical temperature Tc, as a function of the 
parameters of the system, as obtained by using the finite temperature extension of the Slave Boson Mean Field 
Approximation. 

PACS numbers: 73.23.Hk, 72.15.Qm, 73.63.Kv 



I. INTRODUCTION 

The Kondo regime in strongly correlated mesoscopic struc- 
tures has been extensively studied since its observation in a 
single quantum dot (QD) connected to metallic leads ^ The 
increasing interest in investigating the underlying physics of 
these structures is motivated by potential technological appli- 
cations, such as the design of single electron transistor devices 
based on the Coulomb blockade effect, as observed in sin- 
gle QDsj^ A system composed of two tunnel coupled QDs 
connected in series with metallic leads [henceforth called a 
double QD (DQD)] is particularly interesting because it is the 
simplest geometry in which the interplay of two energy scales 
determines its physical properties: (;) the Kondo correlation 
between the spin of each QD and the spins of the conduction 
electrons, and ( ii) the antiferromagnetic correlation between 
the spins of the two QDs. Each of these two regimes prevail 
in different regions of the parameter space and compete in the 
crossover region. As discussed below, these regimes are man- 
ifested very clearly in charge transport measurements. 

Several experimental and theoretical works have appeared 
in the last decade studying DQDs. The continued interest in 
DQDs stems from the early recognition^^ of a possible Non 
Fermi Liquid (NFL) Quantum Critical Point (QCP) separat- 
ing a Fermi Liquid (FL) local singlet antiferromagnetic-phase 
from a Kondo screened FL Kondo-phase in the Two Impu- 
rity Kondo model (TIKM)j^ Subsequent Numerical Renor- 
malization group (NRG)^ calculations on the Two Impurity 
Anderson model (TIAM) detailed the properties of this NFL 
QCP)^ but already pointed out that the inter-impurity hop- 
ping suppresses the critical transition.'' The great flexibility of 
mesoscopic systems, especially semiconducting lateral QDs, 
where a continuous tuning of the relevant physical parame- 
ters is possible, lead to a series of detailed theoretical stud- 
ies of DQDs as a prototype for the TIAM and TIKM. Sev- 
eral papers used different flavors of the Slave Boson formal- 
ism to analyze either the TIAMiiiri^ and the TIKMJ^ De- 



tailed studies have also been done using NRGr«a the Embed- 
ded Cluster Approximation (ECA)ji^ and the Non-Crossing 
Approximation.-" The results obtained confirmed that when 
the even-odd parity symmetry is broken, the critical transi- 
tion is replaced by a crossover. In addition, in the TIAM, it 
was found that, as the inter-dot hopping increases, a coher- 
ent superposition of the many-body Kondo states of each QD 
(forming bonding and antibonding combinations) results in a 
splitting of the Kondo resonance, which leads to a splitting of 
the zero bias anomaly in the differential conductance i ' ^- 
These many-body molecular states should not be confused 
with the single-particle molecular states (separated by an en- 
ergy equal to twice the inter-dot hopping). In reality, the co- 
herent single-particle molecular states were the first to be ob- 
served by R. H. Blick et al:s^ through careful analysis of the 
charging diagram of a DQD. Concurrently, by using a combi- 
nation of charge transport and microwave spectroscopy, Oost- 
erkamp et al}^ probed the formation of single-particle molec- 
ular states (which they called 'covalent bonds') in a DQD by 
varying the tunneling coupling between the QDs, showing the 
possibility of controlling the quantum coherence in single- 
electron devices. Subsequently, Qin et al?^ showed that these 
coherent molecular states are robust even when coupled to 
acoustic phonons created in the system. The first observa- 
tion of a coherent Kondo effect in a DQD was reported by 
Jeong et al.w^ Indeed, the splitting of the Kondo resonance 
(into bonding and anti-bonding many-body states) was then 
clearly observed, as had been theoretically predicted ; "-'^-'^ 
However, no single Kondo peak has been experimentally ob- 
served to date, probably due to the small values of inter-dot 
tunnel coupling required, which leads to a very small overall 
conductance at half-filling. 

In addition, several recent papers have been trying to deter- 
mine under what experimental conditions it would be possible 
to observe manifestations of the NFL QCP mentioned above. 
In the theoretical side, Affleck et al^ have used bosoniza- 
tion to determine the energy scale below which the RG flows 
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to a FL QCP (away from the NFL QCP, as discussed above 
- see also Ref. 1261) . More importantly, an exact functional 
form for the conductance in the crossover region from NFL 
to FL was derivedf2£ opening the doors for its possible exper- 
imental observation. Note that a very interesting experiment 
involving carefully approaching a Cobalt atom in an STM tip 
to another Cobalt atom laying on a gold surface has nonethe- 
less failed to observe the NFL QCP physicsjS Also, M. Lee 
et al. have found through NRG a conduction-band mediated 
superexchange J/, which competes with the direct superex- 
change term Ju ~ At^ /U and dominates for large values of 
U (intra-dot repulsion, and t being the inter-dot hopping ma- 
trix element). These two terms (related to inter-lead charge 
transfer) are known to destabilize the NFL QCP, by breaking 
parity. Indeed, as mentioned above,'' this QCP is closely asso- 
ciated to the TCK fixed point which is believed to have been 
observed in only one system.— In a recent paper, F. W. Jay- 
atilaka et al.r- using NRG, carefully analyze the experimental 
possibility of observing this QCP in a DQD. 

In this work, we study the properties of a DQD (see Fig. [Til 
using the Slave Boson Mean Field Approximation (SBMFA) 
at finite U. The main property we are interested in is the be- 
havior of the so-called Kondo cloud in a DQD. The Kondo 
cloud can be understood as the spatial region occupied by the 
conduction electrons that are collectively involved in screen- 
ing the impurity spini22i^ In the case of the DQD analyzed 
here (Fig. [T]i, at half-filling, when the tunnel coupling be- 
tween them is weak (as compared to the coupling of each 
QD to its adjacent lead), each QD will form its corresponding 
Kondo cloud with the conduction spins of the lead to which 
it is directly connected. However, as the inter-dot coupling 
increases, the system is driven through the crossover region 
between the Kondo regime and the antiferromagnetic 'molec- 
ular' regime. When passing through this region, the Kondo 
cloud associated to each dot should 'shrink' and disappear ac- 
cordingly. 

In a ID system, the Kondo cloud is characterized by its 
length, which can be estimated by assuming that the mean free 
path of the many-body quasi-particle is related to a time scale 
associated to the Kondo temperature Tk- If one assumes that 
the electrons scattered by the single impurity propagate with 
the Fermi velocity vp, the Kondo cloud length, ^, is estimated 
asiii 



where h and are Planck and Boltzmann constants, respec- 
tively. Therefore, in the present case, it is important to deter- 
mine whether for the DQD the Kondo cloud length also scales 
with the inverse of Tk, as shown in Eq.[T] 

From a theoretical point of view, the Kondo cloud has been 
studied through different approaches.ii The most common 
one is based on the study of the dependence of the spin-spin 
correlation with the distance between the impurity and the 
conduction electrons.^- Using a variational approach it has 
been suggested that ^ does not play a significant role in the 
physics of a system of impurities in 2- and 3 -dimensions f22, 
However, the Kondo cloud is an important concept to analyze 



FIG. 1: (Color online) The figure sketches an artificial molecule con- 
sisting of two QDs with intra-dot Coulomb interaction U , and tun- 
nel coupled to each other through the inter-dot matrix element tap. 
In addition each QD is tunnel coupled to its adjacent metallic lead 
through the matrix element The leads L and R act as charge 

reservoirs, being in thermodynamic equilibrium with the DQD. 



the conductance properties of 1 -dimensional systems, like the 
one studied in this paper 

Very recently, the authors studied the behavior of this cor- 
relation for arbitrary distances from a single impurity f^i This 
study has been done by analyzing the effect of the Kondo res- 
onance on the local density of states (LDOS) away from the 
impurity. Using this approach, it was possible to show that ^ 
behaves according to Eq. [T] As a consequence, it is in prin- 
ciple possible to determine Tk by studying the length of the 
associated Kondo cloud. 

To be more precise, the main objective of this paper is to 
answer the question: how does the Kondo temperature and 
the Kondo cloud depend on the parameters of the DQD? To 
this end, we employ the finite-U SBMFA, which requires 
a more involved numerical calculation than the infinite-U 
SBMFAriii^ and therefore eliminates the misconceptions cre- 
ated by the artificial introduction of an extra antiferromagnetic 
inter-dot interaction J in the Hamiltonian. We also study the 
system at finite temperature and estimate Tk by associating it 
to the temperature above which the width of the Kondo peak 
in the LDOS vanishes. In this limit, within this approxima- 
tion, the QDs decouple from the rest of the system. We then 
compare the Kondo temperature obtained from the length of 
the Kondo cloud with the one obtained from the slave-boson 
criterion just described above. The results agree very well 
in the region of parameter space where the system is in the 
Kondo regime. 

The paper is organized as follows: In Sec. II we present the 
model for the DQD and describe the finite-U SBMFA, while 
in Sec. Ill we present SBMFA results at half-filling. Sec. 
IV is dedicated to the study of the Kondo cloud inside the 
leads, presenting a way of calculating its extension and from 
it the Kondo temperature. In Sec. V we compare the results 
obtained for Tk from the calculation of the Kondo cloud ex- 
tension with those obtained by calculating the width of the 
Kondo peak at finite temperature. In Sec VI we present our 
conclusions, and relegate to the Appendix some more techni- 
cal results. 
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II. MODEL AND SLAVE BOSON MEAN FIELD 
APPROACH 

The system of two quantum dots presented in the Fig. [T] is 
described by an Anderson Hamihonian which can be sepa- 
rated in three parts, H = Hq + Ht + HieacU, where 



E 



z— a,/3 



(2) 



describes the isolated QDs, in which {ci„) creates (anni- 
hilates) an electron with energy and spin a in the ith QD 
(i = a, /3), and the second term corresponds to the Coulomb 
interaction U in each QD. 

(T 

+ ^ tap {c'aaCfia + H .c) (3) 
a 

describes the connection of each QD with its adjacent lead 
(first term) and the tunnel coupling between the QDs (second 
term). Finally, 
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H.c. 



(4) 



describes the leads, modeled as two semi-infinity tight- 
binding chains of non-interacting sites connected through the 
hopping term t. For simplicity, we take the applied gate poten- 
tial to be equal in both QDs, i.e., £0 = 6/3 = Vg, and we con- 
sider symmetric coupling to the leads, = tn = t' . Here- 
after, we choose t = 1 as our energy unit and set ft = fc^ = 1 . 

Within the slave boson formalism, the physics underlying 
the Kondo regime is brought into the model by the introduc- 
tion of the auxiliary bosons e^, pi^r, and die, which project the 
Hilbert space onto space sectors with zero, one, and two elec- 
trons, respectively. To accommodate the new bosonic field 
with these operators, the Hilbert space is naturally enlarged, 
and the single electron operator Cjcr is replaced by a 

quasi-electron operator Zcia {Z^c\^), where the Z-operator, 
in mean field approximation becomes just a real factor— 

= [1 - {d,f - {p,,f]~^/\{e,)lj>,„) + X 
x[l-{e,f -{p,,f]-^/\ (5) 

where i denotes the ith QD and {e.i), (pia), and {di^) the mean 
value of the slave boson operators. The full Hilbert space has 
to be restricted to the physically meaningful sector by impos- 
ing the constraints 



and 



^(p,.>' + (d.)'-l = 0, 



riia - ipia)'^ - {dif = 0, 



(6) 



(7) 



via Lagrange multipliers A^-* and A^-^^ Considering the hy- 
bridization of the fermion operators in Hq and Ht, and also in- 
troducing the Lagrange multipliers A^^ and \\^, we can write, 
the effective Hamiltonian as. 



H, 



J2 ^^(^^ 

i—a.,P 



E 

RU=f>) 



trZ 



cLcj^ + H.c. 



E 

i—a.,13 



E^2. 



^tapZ^ [c^aaCfia + H .c] 
a 

e.>' + E^^'->' + <'^')'-l + 

a 



(8) 



where 6^ = 6^ + X^^ is the renormalized quasi-fermion energy. 
It is important to notice that the replacement of the single- 
fermion operator, in mean-field approximation, results to be 
equivalent to a renormalization of the connections t' and t^p 
by a multiplicative parameter Z and Z^, respectively. With 
the effective Hamiltonian Hf.f / we can obtain the free energy 
of the system, which is minimized with respect to each com- 
ponent of the set of parameters. 
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{Ca, e/3,Pa<T,P/3<T, da, dp, A", Aj* , Aj^ , A^^ } . (9) 



The minimization of the free energy provides a set of non- 
linear equations that has to be solved in a self-consistent way. 
See Appendix B for details. 



III. DIFFERENT MANY-BODY REGIMES 

In this section, we present SBMFA numerical results for the 
low-temperature physics of the DQD shown in Fig. L These 
results were obtained in the two different regions of the pa- 
rameter space, as described in the Introduction, and allow the 
analysis of transport properties in the two well characterized 
quantum states of the DQD, as well as in the crossover re- 
gion between them. In addition, some results for the molec- 
ular Kondo regime, which is characterized by an effective 
Coulomb interaction Uejj-, will be presented as an Appendix. 

A. Atomic Kondo Regime at Half-Filling 

To start the discussion concerning the different regimes of 
the system and the crossover between them, we present in 
Figs. 121 El and m the renormalized energy level ia{p) and the 
local density of states (LDOS) of the QDs, the conductance G, 
and the renormalization parameter Z^, respectively. All these 
quantities (except for the LDOS in Fig. |2lB) are calculated 
as a function of the gate potential Vg and for different values 
of tap, and Ep = in both reservoirs. If the gate potential 
satisfies Vg sa —U/2 (i.e., the DQD is occupied by two elec- 
trons), then, for small values of tap, the system is expected to 



4 



be in the Kondo regime resulting from the singlet state cre- 
ated by the antiferromagnetic correlation of each QD with 
the conduction electron spins of the corresponding lead. The 
formation of this Kondo state is characterized by the plateau 
structure, around Vg ~ —U/2, at the Fermi level, observed in 
the results in Fig. I2A) for taf) = 0.025, for example. This 
plateau is related to the opening of a conducting channel (at 
u! = 0), which allow charge transport through the DQD. As- 
sociated with this conduction channel there is a Kondo reso- 
nance at the Fermi level, which can be clearly seen in Fig.|2lB) 
[(black) square curve, for t^p = 0.025], which shows re- 
sults for the QDs' local density of states (LDOS) close to the 
Fermi level for Vg = —U/2. Note that the Kondo peak is 
split (and suppressed) for the smallest value of tai3 = 0.025 
in Fig.|2B), and that this splitting becomes more pronounced 
as tap increases!^ Consequently, the conductance (shown in 
Fig. |3]l is also suppressed from its maximum value G = Go 
atVg = —U/2. As mentioned in section II, another important 
characteristic of the Kondo regime, in the finite-U slave bo- 
son approach, is the renormalization of the couplings through 
the parameter Z (between both QDs and between each QD 
and its adjacent lead), as shown in Fig. |4] Note that for 
ta0 = 0.025 [(black) square curve], this renormalization is 
the strongest. This result shows a large reduction of the pa- 
rameter Z, reaching Z'ft!v'0.05«0.22, strongly suppressing 
the renormalized hoppings t' and tap, by factors Z and Z^, 
respectively [see eq. (8)]. This is an indication of the forma- 
tion of a Kondo state, where spin fluctuations are enhanced 
and charge fluctuations in the QDs are suppressed, which is 
reflect in the decrease of the effective hoppings connecting 
the QDs to each other and to the leads. One should remark 
that the split Kondo peak [(black) square curve] in Fig. |2B) 
is associated to the many -body coherent states which had been 
theoretically predicte d"''^''^ and experimentally observed.— 



B. Molecular Regime at Half-Filling 

With the gradual increase in the magnitude of the connec- 
tion tap, an antiferromagnetic interaction J = tap"^ /U de- 
velops between the spins localized in each QD, competing 
with the Kondo state and, for large values of tap, it is re- 
sponsible for the total suppression of the Kondo state at half- 
filling. From Fig. EtA) we note the progressive destruction, 
as tap increases, of the plateau in the la(i3) vs Vg curves, 
around half-filling. Simultaneously, it is possible to observe 
the development of two other plateaus, for Vg/U > and 
Vg/U < —1, that correspond to the emergence of one- and 
three-electron Kondo regimes, respectively, which will be dis- 
cussed in Appendix A. In Fig.|2B) (where Vg — ~U/2, i.e., 
particle-hole symmetric point), the transition to the molecular 
regime, accompanied by the destruction of the Kondo reso- 
nance [(black) squares curve), and the formation of the molec- 
ular antiferromagnetic state [(blue) triangle curve) is reflected 
in the LDOS Pa(p) of the QDs. Note in Fig. ^B) that for 
tap = 0.125 [(green) down-triangle curve], the separation be- 
tween the peaks is already w 2tap, indicating that these peaks 
are associated to the coherent single-particle molecular states. 
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FIG. 2: (Color online) (A) behavior of the renormalized local energy 
state €a(p) of the QDs as a function of the gate potential Vg. The 
parameters used are U = 0.5, t' = 0.15, Fermi energy Ef = 0, and 
with different magnitudes for tap (see legend). As tap increases, the 
plateau observed at eci(/3) = 0, characteristic of the Kondo regime, is 
gradually suppressed. Note that for the largest value of tap = 0.45 
one has now two plateaus located at Sa(p) ~ itap, which (as noted 
in Appendix A) are associated to the molecular Kondo regimes at 
quarter-filling. (B) QD's LDOS for Vg = ~U/2, and the same 
tap values as in panel (A). For the smallest value of tap ~ 0.025 
[(black) squares curve) a splitting of the Kondo resonance is already 
observed. Note that for tap = 0.125 [(green) down-triangles curve], 
the separation between the peaks is already close to 2tap, which 
is the separation expected for single-particle molecular orbitals, as 
mentioned in the Introduction. 



as experimentally observed in Refs. l2ll - l23l The transition to 
the molecular regime can also be observed in Fig. H] show- 
ing that the renormalization of the dot hopping connections is 
reduced, as the renormaUzation parameter rapidly approaches 
the value Z«1.0 as tap increases. This effect is reflected in 
the conductance of the system as presented in Fig. [3] which is 
strongly suppressed at half-filling for larger values of tap . The 
plateaus outside the half-filling regime shown in Fig. 2(A), for 
the largest values of tap correspond to the conductance of the 
DQD at the molecular one- and three-Kondo regime that, as 
mentioned, are discussed in Appendix A. 
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FIG. 3: (Color online) Conductance as a function of Vg for t^p — 
0.025 [(black) squares curve], 4^/3 = 0.05 [(red) circles], tafi — 
0.075 [(blue) up- triangles], ta/s = 0.125 [(green) down- triangles], 
tai3 = 0.25 [(magenta) diamonds] and t^is = 0.45 [(light green) 
left- triangle]. The other parameters are t' = 0.15 and U — 0.5. 
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current work we consider a similar expression. 



FiN) 



pf{uj)]p^i^p){uj)du, (10) 



to quantify this perturbation. Instead of using a Gaussian 
function with width Tk to convolute the density of states and 
eliminate the Friedel oscillations?^ allowing the study of the 
region near the Fermi level, we use the LDOS Paip) of the 
QD. Note that, as pa = pp, we will omit the subscripts a 
and P from the rest of the equations in this subsection. The 
use of this definition for the convolution function (p) has the 
advantage of directly incorporating into F{N) the physical 
information associated with the Kondo ground state. In this 
expression, p^ and p"^ represent the LDOS calculated in the 
A^th site inside the lead (L or R) with the system in and away 
from the Kondo regime, respectively. The last condition is 
enforced by disconnecting the DQD from the leads, i.e., by 
calculating the LDOS for t' = 0. It is well known that, within 
the SBMFA, this is equivalent to taking the system to a tem- 
perature T > Tk-^ 

The LDOS appearing in the integrand of F{N) is propor- 
tional to the imaginary part of the Green's function Gnn de- 
fined at the A^*'' site inside the semi-infinite metallic leads. 
Considering the left lead, for instance, we can write 



GnN — 9NN + 9Nlt'GaN- 



(11) 



In this expression gj^N is defined as the Green's function at 
site N of the left lead when t' = and g^i satisfies the equa- 
tion 



9N1 - t ffL : 



where 
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(12) 
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FIG. 4: (Color online) Renormalization parameter as a function 
of Vg for the same t^^, U, and t' values as in Fig.[3] The inset shows 
the small renormalization of Z characteristic of the molecular Kondo have 
regime. 



is the Green's function defined for the left lead. For G 



GaN — Gnc = gNit'Go 



aN{Na} 



(14) 



IV. KONDO CLOUD 



where Gaa is the dressed function at the QD a. Substituting 
GaN into eq. (11), yields 



A. Cloud Extension Function 



GnN — 9NN + t''^{gNlYGa 



(15) 



In this section we introduce a new perspective to the Kondo 
problem in a strongly coupled DQD at half-filling, based on 
the analysis of the extension of the Kondo cloud inside the 
metallic leads. To analyze the extension of the Kondo cloud 
we use the method developed in Ref. [sil where the authors 
analyze the propagation into the leads (away from the QD) of 
the perturbation in the LDOS, introduced by the Kondo reso- 
nance at the QD. To study this propagation, a function F{N) 
was defined, which quantifies the perturbation produced by 
the presence of the Kondo cloud in the iVth site of the lead 
(counted from the border of the semi-infinity chain). In the 



We note that the Kondo physics is introduced into Gnn 
through the term proportional to t'^. The imaginary part of 
this function is proportional to the LDOS p\j at the iV*'' site. 
The non-Kondo solution is obtained by eliminating the ef- 
fects resulting from the presence of the QD by considering 
t' = in the expression of Gnn- So, p^^ is defined as 
p'^*-^ = -(l/7r)lm5ArAr. 

Considering the analytical expression for and p^^' we 
obtain 



fe _ nk 

Pn Pn 



1 



Im[t'2(gjvi)'G„J. 



(16) 
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Substituting gpfi into eq. ( fTSI l we obtain 

„nk ^ T„\+I2~2N 4-2N-2r' 1 t\n\ 

Pn - Pn = ^ CraaJ- (17) 

TT 

This expression, when substituted into eq. ( fTOl i. results in 

2 /■ + 00 



F(iV) 



Imrgri'''"'G„„(^)]p(^)'^w, (18) 



where all the effects of strong correlation present at site N are 
contained in the Green's function of the QDs. Note that when 
tap = one should obtain the results from previous calcula- 
tions for a single-QD Kondo cloud.- For finite tap, the com- 
petition between Kondo and antiferromagnetism is contained 
in Gaa{^), and therefore should be reflected in F{N). 



B. LDOS 

The LDOS p% and p'ff^, calculated at tiie site = 50 of 
the leads, for the system with tap = 0, are shown in Fig. |5] 
As in the rest of this paper, the results presented in this fig- 
ure were obtained for Vg = -U/2, U = 0.5, and t' = 0.2. 
The LDOS p^v^ was calculated for an isolated lead {t' ~ 0). 
The Kondo resonance at the QD is shown in the (blue) tri- 
angle curve of Fig. |5] On the other hand, the LDOS /0^^5q 
[(red) circle curve] shows a small peak at the Fermi level. The 
presence of this peak is a sign of the Kondo resonance 'prop- 
agating' through the sites of the leads and reflects the exis- 
tence of the Kondo cloud, with an extension ^, dependent on 
the Kondo temperature Tk- In our context, the extension of 
the Kondo cloud is obtained from F{N) (see next section for 
more details). It is important to note here that the propagated 
Kondo peak along the leads appears as a resonance or an an- 
tiresonance depending on whether the site N is even or odd, 
respectively!^ 

Figure |6] shows the behavior of the LDOS for different 
values of tap, illustrating how the competition between the 
Kondo regime and the antiferromagnetic correlation manifests 
itself in the LDOS at an arbitrary site {N = 50). The (ma- 
genta) dashed curve is the LDOS for a lead with no connection 
to the dots, which, as a reference, corresponds to the expected 
result obtained for a non-Kondo regime. The (black) down- 
triangle curve corresponds to tap — {f ~ 0.2) and shows 
the characteristic Kondo resonance propagated to = 50. It 
is clear from the figure that an increase of the tunnel coupling 
tap drives the system from a Kondo regime to a non-Kondo 
ground state, as the LDOS at (and around) the Fermi level 
is clearly suppressed, approaching the value for the discon- 
nected DQD [(magenta) dashed curve]. 



V. KONDO TEMPERATURE 

A. Kondo Temperature from the Cloud Extension Function 

Figure|2]shows the logarithmic dependence of the function 
F{N) for various inter-dot couplings tap- Differently from 




FIG. 5: (Color online) The figure shows the effect produced by the 
presence of the impurity a(/?) in the LDOS calculated at the site 
A'' = 50 [inside the lead L(R)]. The (black) squares curve shows the 
LDOS (at site A'^ = 50) for the isolated lead; the (red) circles curve 
shows the LDOS when each lead is connected to its respective QD, 
for tap = 0; and the (blue) triangles curve shows the LDOS for the 
impurity. 
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CO 

FIG. 6: (Color online) The figure shows the effect produced by the 
antiferromagnetic interaction J — At^p/U between the QDs in the 
LDOS at site AT = 50 inside the leads. The (black) down-triangle 
curve shows the LDOS calculated at A'' = 50 for tap = 0. The (red) 
up-triangles, (blue) circles, and (cyan) down-triangles curves are ob- 
tained for tap ~ 0.03, tap ~ 0.06, and tap = 0.09, respectively. 
The (magenta) dashed curve shows the LDOS for the isolated leads 
{t' = 0). The other parameters are Vg = ~U/2, U — 0.5, and 
t' = 0.2. 



the result obtained for a system of one QD connected to a 
metallic lead;^ this function presents now an oscillatory be- 
havior with a frequency associated to the inter-dot coupling. 
However, similarly to the single-QD case, the extension ^ of 
the Kondo cloud can still be obtained, in this case by the expo- 
nential decay of its envelop function. The physical informa- 
tion contained in this function is extracted from the straight 
lines, tangent to the logarithm of the F{N) function. Specif- 
ically, the extension of the Kondo cloud, as a function of the 
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FIG. 7: (Color online) Natural logarithm of the cloud extension func- 
tion F vs the distance N (from the border of the semi-infinite chain) 
for ta/s — 0.0 [(black) squares], ta^ — 0.04 [(blue) triangles], and 
tai3 ~ 0.07 [(red) circles]. The straight lines are tangents to InF and 
correspond to tap = 0.01 [(purple) line], ta0 = 0.02 [(black) line], 
tap = 0.03 [(red) line], tap ~ 0.12 [(green) line]. The other 
parameters are Vg = -(7/2, U = 0.5, and t' = 0.2. 



inter-dot connection, can be obtained from the slopes of these 
lines, which are proportional to 1/^ (see eq. (23) in Ref.lsH). 

In Fig. |2] we show the function In F{N) for tap = 
[(black) squares curve], tap = 0.04 [(blue) triangles curve], 
and tap = 0.07 [(red) circles curve], as well as their respec- 
tive tangent lines. For intermediate values of tap we present 
only the tangents to In F{N). As mentioned above, the slopes 
of these straight lines are proportional to 1 and allow us to 
obtain the Kondo temperature through the expression ^ = 
proposed in Ref. |3ll The values obtained for the Kondo tem- 
perature Tk are presented in the (red) circles curve in Fig.|9] 
as function of tap (we label it as T^'~^, i.e., the Kondo tem- 
perature obtained through the extension of the Kondo cloud). 

We note that, for —ff^ < Tk, which implies tap < 0.07, Tk 
presents an exponential behavior, in accordance with results 
obtained by Aono and Eto^ using a slave boson formalism 
with an infinite Hubbard U. 




of temperature T for values of tap ranging from 0.00 [(black) left- 
most curve] to 0.13 [(purple) rightmost curve]. The other parameters 
are t' = 0.2, U = 0.5, Vg = -U/2, and Ef = 0.0. 



from this result that the parameter vanishes rapidly when 
the temperature approaches the characteristic value Tc- 

The values of obtained for the different values of tap are 
represented by the square (black) curve in the Fig. |9](we la- 
bel it T^"*), and it agrees with the values obtained for T^^ 
from the function \nF{N) associated to the extension of the 
Kondo cloud in the metallic leads. We observe a qualita- 
tive and semiquantitative agreement between these two ways 
of calculating the Kondo temperature, giving support to the 
interpretations we are proposing. We note that increasing 
tap > 0.07 the values of Tc become increasingly different 
from the Tk values obtained from the extension of the Kondo 
cloud. These are precisely values for which tap'^/U > Tk, 
and the discrepancy starts at tap ^ 0.07, where the DQD 
enters a crossover region which extends up to values of tap 
for which 4:tlp/U > 2T^\ where the DQD enters into an 

antiferromagnetic phase (T^^ is the Kondo temperature of 
each op when tap = 0.0), according to previous works (see 
Ref. ITU and references therein). 



VI. CONCLUSIONS 



B. Kondo Temperature in the Finite Temperature SBMFA 

In this section we determine the Kondo temperature Tk 
of the DQD by extending the slave boson formalism for fi- 
nite temperature. In this formalism, the impurity is decou- 
pled from the rest of the system for a critical temperature Tc. 
Although phase transitions are typical artifacts of mean field 
solutions, the temperature for which this decoupling occurs 
can be considered as a reasonable approximation to the Kondo 
temperature, which physically corresponds to a crossover be- 
tween two regimes. The decoupling occurs through the pa- 
rameter Z, which renormalizes the coupling of the QDs to the 
leads, t' = Zt' and tap = Z'^tap, which is vanishingly small 
when T w Tk- The parameter Z, as a function of temper- 
ature, is studied in Fig. [8] for various values of tap- We see 



In this paper, the interplay between the antiferromagnetic 
interaction and the Kondo effect for a DQD at half-filling has 
been studied in detail. In particular, we have analyzed the 
dependence of the extension of the Kondo cloud on the inter- 
dot tunnel coupling (tap) when the system is driven from the 
half-filling Kondo regime to the antiferromagnetic molecular 
ground state, as tap increases. The extension of the Kondo 
cloud and its associated Tk were obtained by analyzing the 
propagation inside the leads of the Kondo resonance 'origi- 
nating' at the dots. In the Kondo regime, the Tk results ob- 
tained through the Kondo cloud extension were almost iden- 
tical to the 'Kondo transition temperature' Tc obtained from 
the finite-temperature extension of the MFSBA. Although we 
know that the decoupling of the DQD from the leads (occur- 
ring at Tc) is an artifact of the MFSBA, our results confirm 
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FIG. 9: (Color online) Kondo temperature T^"^ and Tj?"* as a func- 
tion of ta0. The (red) circles curve corresponds to Tk obtained from 
the extension ^ of the Kondo cloud, while the (black) squares curve 
corresponds to Tk as obtained through the critical temperature Tc 
for which the QD's are decoupled (see text), with Z ^ 0. The insets 
show the exponential fits obtained for each Tk in the region where 
the system is expected to be in the Kondo regime. The values for the 
other parameters are Vg = ~U/2, U — 0.5, and t' — 0.2. 



that the Tc can be reliably taken to be « Tk in a DQD, very 
much in the same way as in the case of one impurity. Note 
that the dependence of Tk obtained by both approaches agrees 
with previous results by Aono and Eto. - 

The study of the Kondo regime in the DQD was done within 
the context of a finite U treatment. As U is finite, the Hamil- 
tonian incorporates the correlation between the QD's spins. 
This allows us to avoid introducing an artificial antiferromag- 
netic interaction between the QD's, normally incorporated in 
U = oo approaches (as done, for example, in Ref. 16), which 
provides more confidence to our results. Its important to note 
that the study developed in this paper opens a set of concep- 
tual ideas that can be useful to understand the effects of the 
RKKY interaction between two QD's located at an arbitrary 
distance between themselves, as well as the transport associ- 
ated to them. 



Appendix A: Molecular Kondo Regime at Quarter-Filling 

Given the importance in a DQD, as mentioned in the In- 
troduction, of the concept of molecular states (or molecular 
orbitals), the authors will briefly describe the properties of the 
molecular orbitals in a DQD, as done previously for related 
systems/^ For recent work by some of the authors, see Ref- 
erences l36ll37[ A system of two or more coupled QDs is said 
to be in the molecular regime when the charge transport oc- 
curs through orbitals that involve linear combinations of lev- 
els in each QD (which could be called 'atomic' orbitals). The 
molecular Kondo effect occurs when an unpaired spin residing 
in a molecular orbital is screened by conduction electronsi^ It 
is clear that, in the molecular Kondo regime, the charge at each 
QD in a multi-QD system does not correspond to an integer 



FIG. 10: (Color online) Conductance as a function of tafj for 
Vg/U = -0.5 [(black) squares], Vg/U = -0.14 [(red) circles], 
Vg/U = 0.0 [(blue) up- triangles], Vg/U = 0.08 [(cyan) down- 
triangle], Vg/U = 0.22 [(magenta) diamond], Vg/U = 0.6 [left- 
triangle (gold)] and Vg/U — 1.2 [(royal) right- triangle]. The other 
parameters are t' = 0.15 and U = 0.5. 



value. As a consequence, the non-integer QD charge ground 
state of a multi-QD structure cannot be used as a criterion to 
identify the system as being in a fluctuating valence regime, 
as it could be in a molecular Kondo ground state, involving 
more than one QD. The determination of the real nature of 
the ground state of these systems requires a careful analysis, 
mainly for parameter values where the system is in a transition 
between the molecular and the atomic regimes. 

In the case here analyzed for a DQD, by adjusting Vg to 
have one (or three) total number of electrons in the DQD, 
and considering relatively large values of tajs, the quarter- 
filled molecular Kondo regime is accessed. In the context of 
the finite-?/ slave boson formalism, this regime is character- 
ized by two plateau structures in £a{i3) [as observed at half- 
filling in Fig.|2jA)] at ea(^) = itai3- The formation of these 
plateaus is associated with two new resonances which allow 
charge transport through the DQD, as can be seen by the well 
separated double peaks in the conductance in Fig. [3] [(yellow) 
left-triangle and (magenta) diamond curves, for tap = 0.45 
and 0.25, respectively]. 

Figure|4]shows that as the system is driven out of the Kondo 
regime, at Vg = ~U/2, by increasing tap, the parameter 
Z — > 1.0, thus eliminating the renormalization of the hop- 
ping matrix elements [see eq. (8)]. As an important character- 
istic of this result, we note two small suppressions of for 
the (yellow) left-triangle and (magenta) diamond curves in the 
one- and three-electron regions of Vg. This region is ampli- 
fied in the inset of Fig. |4] and reflects a small renormalization 
provided by the Z parameter, characteristic of the molecular 
Kondo regime. 

To investigate in more detail the effect of the different 
regimes in the transport properties of the DQD, in Fig.fTOlwe 
present the conductance as a function of tap for different val- 
ues of Vg. These results are fitted by the analytic expression 
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Appendix B: The minimization of the free energy in the SBMFA 



fitap) 



i^' + t'lp 



_4g2p2' 



(Al) 



where F = ■Ktj^^^j^-^p{EF) is the coupHng constant and i = 
a{f3). This expression is an extension of the one obtained by 
Georges and Meir— and is valid also away from the particle- 
hole symmetric point, Vg ~ —U/2. Exactly at Vg = —U/2, 
the DQD has one electron in each QD and the conductance 
is maximum when tap = F. For higher values of tap, as 
Vg moves away from — [/ /2, the DQD approaches the quarter 
filling molecular Kondo regime. This occurs when the charge 
of the system is nearly one [Ua + np^l) or three electrons 
(ua + np^3). This regime in the conductance is reflected in 
Fig. |3]by the two peaks which occur near the region of one 
and three electron occupation. Figure [TO] shows the conduc- 
tance as function of tap for various values of Vg, that corre- 
spond to the positions of the conductance peaks in Fig. |3] (for 
Vg> — U/2). From the results obtained for > one can 
identify two crossover regions for each Vg as tap increases: 
one corresponding to the transition from an empty dot situa- 
tion to a molecular Kondo regime, with an increasing conduc- 
tance region and the other, for values of tap above the maxi- 
mum of the conductance, corresponding to the antiferromag- 
netic molecular regime where the conductance decreases, and 
the DQD's occupancy increases. 

To better understand the physics underlying the strong 
coupled QDs in the transition region between the different 
regimes in the parameter space defined by tap and Vg, we 
use the exact solution obtained in Appendix D, when the QDs 



are disconnected from their leads. We obtain that E, 



(le) _ 



(2e) 



2en 



eo — tap and Eq'^''' = 2eo H j-j^g jowest en- 
ergies in the one- and two-electron sector, respectively. Con- 
sidering these expressions in a situation in which the charge 
occupation is increasing, we conclude that the inclusion of an 
extra electron into the DQD requires an extra energy given by. 



E, 



(2e) 



2E, 



(le) 



2tap 



— 2 ■ This extra energy 



can be identified with an effective Coulomb interaction 



U^JU 



2t 



aP 



^Gtlp 



which, for small values of U /tap, is w U/2 + {U'^/16tap)- 
The addition of a second electron in the DQD is achieved by 
adjusting the gate potential Vg. In fact, for Vg w tap — U /2, 
the DQD is basically double occupied, characterizing a tran- 
sition from a molecular Kondo regime (at quarter-filling) to 
an antiferromagnetic state, or to a two-impurity half-filling 
Kondo regime, depending upon the ratio / J being less or 
greater than unity, respectively. 



The free energy of the DQD is given by 



-PEihi) 



(Bl) 



where /? = l/ksT, and 7; denotes the Lagrange multipliers. 
Differentiating F with respect to 7; we obtain 



dF 



l3Tki 



d'yi 
d'yi 



(B2) 



where 7; denotes all the components of the 7-set defined in 
eq. (|9l). In addition, we have adopted an approximation as- 
suming that the mean value of the derivative of the internal 
energy with respect to the components ji is approximately 
equivalent to the derivative of the mean value of this energy 
with respect to these componentsi^ 



Appendix C: Conductance calculation 

The conductance is obtained from the Green's functions of 
the system. For the QD's a{(3) we have the local Green's 
functions 



MP0) 



9a{P)c 



l-ap9acrgpa 



(CI) 



obtained through a diagrammatic expansion which incorpo- 
rates in the QD a{/3) the physics underlying the complete 
system, including the electron reservoirs. The function ga{p)a 
that appears in this expression describes the sub-system com- 
posed by the QD «(/?) connected to the L{R) reservoir. The 
calculation of this function results in 



9a(P)o 



9a{P)o 



(C2) 



1 - t''^9a(P)cr9L{R)cr ' 

where ga{p)a is the single-particle Green's function associ- 

(A2) ^"^^^ "('^)' ^^^^^ 9L(R)a- = is the 

L(R) reservoir's Green's function projected onto its nearest 
QDa(/3). 

The non-local Green's functions of the system are given by 



G 



aPiPa) 



and 



G 



La(Rp) 



9a{P)atapgp(a)a 
- tip9a{P)agp{a)'y ' 



gL(B.)<jt'ga(p)cr 
1 ~ tlpgaagpa 



(C3) 



(C4) 
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FIG. 11: (Color online) Schematic view of the decoupled DQD pre- 
sented in Fig. 1. 



which are associated to charge transport between the QDs and 
from the L{R) reservoir to the QD a(/3), respectively. This 
propagators can be used to obtain the conductance G of the 
system. Using GJ^^, for example, we obtain from the Keldysh 
formalism^ the following expression for the conductance: 



G 



2c 2 2 
AITT 



-t-oo 



(C5) 



where e is the electron charge, h is Planck's constant, and 
fL{R){^) the Fermi distribution function associated to the 
L(R) reservoir 



where and Ui are respectively the energy of the local state 
and the Coulomb interaction in the i-th QD, tap is the hopping 
term and a the spin projections of the electrons in the QDs. 
The DQD may have an electron occupancy of = 1, 2, 3, 4. 
Thus, considering the system with two electrons, = 2, we 
define the basis 



in,o) 
|o,n> 

^[it,;> 
1 



i;,t)] 



V2 

lt,t) 

i;,;> 



[it,;> + i;,t>] 



(D2) 
(D3) 

(D4) 

(D5) 

(D6) 
(D7) 



where \ipi), \(p2) and lips) are states with total spin 5*7- = 0, 
while \ip4), 1 1^95) and \(Pq) are states with 5^ = 1. Written 
in this basis, the Hamiltonian can be separated into blocks, 
which correspond to the projections 5*^ = and 5t = 1 of 
the total spin. The block-matrix that corresponds to a spin 
projection 5't = 1 is already diagonal and its eigenvalues are 
E4 = = = 2eo (for simplicity we consider ea = = 
eo). These energies are associated to states with 5t = 1 and 
5, = 1,0,-1. 

The block-matrix corresponding to 5*7- = is given by 



Appendix D: Exact Solution of isolated DQD at half-filling 

In order to better understand the characteristics of the an- 
tiferromagnetic and the Kondo molecular regime accessed 
when the DQD is at half- or quarter-filling, respectively, we 
calculate the exact solution when the DQD is decoupled from 
the metallic leads, as shown in Fig. (fTTI) . The exact solu- 
tion obtained in this Appendix allows the identification of 
the regions in the parameter space corresponding to different 
regimes and the crossovers between them. 

The Hamiltonian for the decoupled DQD in Fig. (fTTT i is 
given by 

>=a-^ i=a,l3 

+ W(cLc,3<T +4<TCa<T), (Dl) 



/ 2eo + C/ ~V2tap \ 
Hf^e) = 2eo + [/ V2tap , 

\ -V^tap V^tap 2eo / 

with the following eigenvalues: 

El = 2eo + U (D8) 

E2 = ^[4eo + U+^m + ml^] (D9) 

E3 = ^[4eo + U-^U^ + mtl^]. (DIO) 

These results allow (as shown in the Appendix A) a better un- 
derstanding of the processes underlying the transition between 
the different regimes of the DQD. 
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